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Abstract 

Isotopic spin dependent lattice gas model is used to examine if it produces the isoscaling be- 
haviour seen in intermediate energy heavy ion collisions. Qualitative features are reproduced but 
quantitative agreement with experiments is lacking. 
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I. INTRODUCTION 



A very interesting feature of heavy ion collisions at intermediate energy is that provided 
the experimental set ups are identical,the ratio of isotope yields from two different reactions, 



1 and 2, R2i{N, Z) = Y2{N, Z)/Yi{N, Z) exhibits an exponentia! 



[ii,y,y,Q,3: 



relationship function 



of the isotope neutron number and proton number Z 

R2i(N, Z) = Y2{N, Z)/Y^{N, Z) = C exp{aN + PZ). (1) 

This is called isoscaling and a and (3 are called isoscaling parameters. Various theoretical 
models have been examined for at least an approximate emergence of this law. A grand 
canonical model for multifragmentation naturally leads to this law but it is also seen to 
emerge as an approximate formula from a canonical model 6l, an expanding excited source 
model [21] and even anti-symmetrised molecular dynamics calculations [7|. In this note we 
investigate how well isoscaling is obeyed in the lattice gas model (LGM). Part of what we 
do has already been looked at in a recent publication 8| but both our perspectives and the 
lessons we derive from our study are different. Although highly schematic, LGM has been 
profitably utilised for investigating phase transition properties and many features of inclusive 
cross-sections. One might expect that isoscaling is a general feature, not depending upon 
finer details of models and hence should be seen in the LGM multifragmentation model. 
The advavantage of LGM is (1) it takes into account the interaction between nucleons, both 
Coulomb and nuclear exactly (although the nuclear interaction is very schematic); (2) it 
takes into account excluded volume effects exactly and (3) decay of hot nuclei need not 
be considered as the cluster algorithm recognises only particle stable clusters. Many other 
models can not handle the above three issues easily although there are other many virtues 
in these other models. 

The methods for calculating fragments in an isotopic spin dependent LGM have been 
described in many places (see, for example, [o], |lo|, 11, 12 1). The nearest neighbor bond 
between unlike particles e„p is set at -5.33 MeV (to produce binding energy of 16 MeV per 
particle); and the bonds between like particles e^p and e„„ are set at 0. The justifications 
for these choices are explained in Coulomb energy between protons is included. 

These specifications are enough to allow one to calculate populations of all fragments 
given the lattice size, the mass and charge numbers of the dissociating system and the 
temperature. But for later purposes we need to estimate the value of symmetry energy 



implied in the model. The reason for this is the following. The parameter a (and also /3 
of eq.(l)) depends upon the value of NjZ of the the dissociating systems in reactions 1 
and 2. Moreover, it also depends upon the value of symmetry energy which adds a term 
Cs(A^ — Z^' I A to binding energy in the liquid-drop formula. Indeed it is a standard practice 
to use a relation ^] between a and the symmetry energy coefficient c^, as: 

« = ^[(^M)?-(^M)I- (2) 

It is therefore imperative to estimate a value of implied in the LGM. 

II. ESTIMATION OF IN LGM 

We obtain ground state energies of a large number of nuclei by Monte-Carlo sampling 
at zero temperature. For a given nucleus with mass number A and charge Z this ground 
state energy is denoted by BE{LGM). Let the value of the Coulomb energy which is also 
available from this calculation be E^.- Since by choice the volume energy is —16^4 MeV 
we can try to deduce the value of surface energy coefficient and the symmetry energy 
coefficient by setting 

BE{LGM) -E^ + 16A = a.A^/^ + c,(A^ - Z^/A (3) 

Ideally the values of and Cg should be the same for all nuclei. But because there is no 
good reason why BE{LGM) should obey this parametrisation exactly, values of and Cs 
deduced from the above relation will change from nucleus to nucleus. We can now try to get 
the "best" values by minimising the sum of the squares of the deviation. We chose isotopes 
of some arbitrarily chosen Z's. The fit with the best values of as and Cg is shown in Fig.l. 

The fit in Fig.l appears to be very good but some words of caution are needed. Of the 
two constants we are trying to get, is by far the most important one; a^A^/^ dominates 
over Cs{N ~ Zy / A which is a smaller perturbation. But it is Cs{N — Z)"^ /A which is presumed 
to be the deciding factor for isoscaling. What is left after subtracting the surface term can 
not be fitted by the parametrisation Cs{N — Z)'^/A very accurately. This is demonstrated in 
Fig.2. 

What is reported here is similar to but not identical with the extraction of and Cg in 
jol. However there the symmetry energy derived had a volume part (like here) but also a 
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surface part. We have absorbed here all the effects of symmetry energy using a volume term 
only in order to test how well LGM calculations follow eq.(2) which is based on a volume 
symmetry energy. 



III. COMPARISON WITH SOME DATA 

We will compare our calculations with two sets of data. We first consider ^^^Sn+^^^Sn 
(reaction 1) and ^^^Sn+^^^Sn (reaction 2) central collision data. Experimental data are given 
in Fig.l of [sl. Isoscaling is seen to be well satisfied with a value of a= 0.361. In Fig.3. we 
show calculated results for i?2i where this is the ratio of {n2{N, Z)) and {ni{N, Z)); {tln^z) is 
the average multiplicity of the composite with neutrons and Z protons. The dissociating 
systems are taken to be yl=168, Z=75 for reaction 1 and ^4=186, Z=75 for reaction 2. 
These are the recommended values [3i] after allowing for losses like pre-equilibrium emissions 
etc.. The average multiplicity is calculated from 100,000 Monte-Carlo events. We try 10,000 
switches between two events. Metropolis algorithm is used. The slopes of the ratios of the 
average multiplicities should correspond to to the measured value of the slope of experimental 
Y2{N,Z)/Y^{N,Z) (Fig.l of [3]). 

No basic calculation with the grand canonical model which computes the value of a has 



been reported. Canonical model calculations are quite successful [13|]. SMM calculations 
before decay of hot primaries show isoscaling quite well with a=0.46; but after decay of 
primaries isoscaling is not obeyed to the same precision and the approximate a value changes 
slightly to 0.44 \ 

Results from LGM with different lattice sizes (N) and at different temperatures (T) are 
shown in Fig.3. Isoscaling is obeyed very well though it is not as good as in experimental 
data. In LGM there is no correction due to secondary decay. The clusters calculated are all 
particle stable. Even though isoscaling is only approximately obeyed we deduce an average 
value of a. It is about 0.20 for 8^ lattice and temperature 5 MeV (compared to 0.36 in 
experiment). A notable feature of Fig.3 is that for Z=l, N can be as high as 6 (Fig.3 shows 
results upto 5); Z=l and A^=6 as a stable composite happens with the proton in the central 
cube and six neutrons at the six faces. The binding energy per particle for this "nucleus" is 
5.33(6/7). Of course in the real world such a nucleus does not exist. 

We now examine how well eq.(2) is obeyed in LGM. We can write eq.(2) as a = 
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f{cs,T)g{l,2) where (7(1,2) is just a property of the two dissociating systems and / in- 
cludes all the effects of symmetry energy. Examining the validity of a can be ambiguous 
in our case as isoscaling is not equally good for different Z's. Fig. 3 shows that Z=2 obeys 
isoscaling quite well. Let us limit ourselves to a value for Z=2. We consider T=4 MeV, 
lattice size 8^, keep Zi fixed at 75, Ai fixed at 168. For Z2 fixed at 75 we vary A2 from 
168 to 186 and calculate a as A2 is varied. As a function of (^2/^2)^ the value of a is 
quite linear (Fig. 4) suggesting that the functional form of (7(1,2) of eq.(2) is accurate for 
LGM. The same can not be said about f{cs,T). If we use eq.(2), Zi = Z2 = 75 and 
Ai = 168, A2 = 186 and Cs = 23.4 MeV as obtained from the least square fit in section II, 
the predicted value of a would be 0.86 (as contrasted with a ~ 0.22 as actually given by 
the LGM calculation). A different point of view is sometimes taken. One takes the value of 
a and deduces the value of Cs using eq.(2). In our case from the value of a = 0.22 we would 
then be lead to believe that the value of is way lower, about 6 MeV. We see no reason 
why a value of Cg would so drastically change from about 23 MeV for isolated nuclei to this 
low value at 4/5 MeV temperature in an expanded volume in co-existence with other hot 
nuclei. It is specially hard to understand this in LGM. The composites in LGM can not 
be squished or expanded. At most, some of the composites may be in excited states with 
slightly different number of bonds. But such widely different value for symmetry energy 
appears unlikely. We remind the reader that there are two parts to the calculations. One 
is: given the lattice size, the number of neutrons and protons, the temperature and the 
bond strengths e„p, e^p and e„n, calculation of thermodynamic properties and many particle 
correlations at all levels. Monte-Carlo simulations solve this many-body problem of LGM 
correctly though numerically. Next comes the question: how does one define clusters, given 
that the physics upto this point has been done correctly. We follow the prescription, first 
formulated for LGM in ^], subsequently reformulated with the same result in JJJ, shown 
to be closely equivalent to the one derived in jisl and is now universally used. A practical, 
reasonable but different prescription is not known. Given that the choice of the values of 
Cnp, Cpp and e„„ is very restrictive, we have no freedom to alter anything. We find it easier 
to believe that the function f{cs,T) as implied by eq.(2) is not correct for LGM. 

The other set of data we use is for ^^Ni and ^^Ni on ^Be. We assume that the much larger 
nucleus Ni engulfs Be so that for reaction 1 we have Ai = 67, Zi = 32 and for reaction 2 
we have A2 = 73, Z2 = 32. Experimental data can be found in 16|, ll7|]. The subset of data 
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we use here can also be found in [6| (Figs. 6 and 7). Experimental results for Ni on Be show 
a larger deviation from isoscaling compared to what is seen experimentally for Sn on Sn. 
Fig. 5 shows the results of our calculation. There are significant deviations from isoscaling. 
The calculated value of a hovers around 0.22 (these averages are always ambiguous since 
isoscaling deviations are more significant here) whereas experimental values are around 0.6. 
Thus, as in the case of Sn+Sn, theory underestimates the value of a. 

In the case of Ni on Be, experimental values of y2(^, 2')/Y'i(A^, Z) are available for a 
large range of Z from 1 to 28. Low and moderate values of Z have an effective a which is 
much smaller than those belonging to Z ^ 22. (Data can be found in 6|, llGj.) Such details 
do come out in LGM as well. We demonstrate this in Fig. 6. This is highly satisfying. Grand 
canonical model can not explain this difference although the canonical model model can and 
indeed fits the data very well [g]. 



IV. SUMMARY 



Our first objective was to see if isoscaling is obtainable in LGM. It does appear that 
approximate isoscaling is obtained in the model. We confronted the calculations with two 
sets of experimental data. There are no free parameters in the model. The experimental 
values of a are larger than what the model predicts, by about a factor of 2. But the model 
does reproduce some remarkable features. It did show a significant increase in the value of a 
going from low Z isotopes to high Z isotopes as seen in measurements. It also gave linearity 
of a with A{Z/AY. In view of the schematic nature of the model and obvious drawbacks, 
these successes are quite pleasing. Underestimation of the value of a remains a problem. We 
have checked that reasonable variations of lattice size or temperature will not correct this 
problem. Moderate variations in the values of the bond strengths do not provide enough 
corrections. 
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FIG. 1: For selected isotopes we show plots of BE{LGM) — Ec + 16^ (stars) and compare these 
with asA?^^ + Cs{N — Z)"^ /A where ag = 17.0 and Cg = 23.4. These values are chosen from least 
squares fit. All energies are in MeV. 
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FIG. 2: To test the accuracy of the parametrisation Cs(N — jA we plot BE{LGM) — Ec + 
WA - UA^/^ and compare that with 23.4 * (TV" — Z)'^/A. All energies are in MeV. 
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FIG. 3: These figures are calculated to compare with central collision data for Sn on Sn reported 
in Fig.l of [31]. For reaction 2, the dissociating system is set at A = 186, Z = 75. For reaction 1 
we take A = 168, Z = 75. Isoscaling is approximately obeyed but the deviations are not negligible. 
The "average" value of q is ~ 0.2. Experimantally for this case isoscaling is better obeyed and the 
value of a is s::: 0.34. 
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FIG. 4: For Zi = 75 and Ai = 168 we vary A2 from 168 to 186 for fixed ^2=75. The dots 
are from LGM calculation where a is calculated for Z = 2 as for this isotope isoscaling is well 
obeyed (sec Fig. 3)). The best linear fit is shown. This comparison tests the accuracy of the 
[(Zi/^i)2 - {Z2/A2f] factor of eq.(2). 
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FIG. 5: For reaction 2, the dissociating system is set at A = 73, Z = 32. For reaction 1, 
dissociating system is set at A = 67, Z = 32. Experimental data for Ni on Be can be found in 
and also in \&\, Figs 6 and 7. 
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FIG. 6: Same case as in Fig.5, but we also plot the ratios for much higher Z isotopes. Note the 
much higher value of the slopes for large Z nuclei. Experimental data confirm this. The high Z 
cases are not shown in the lowest right panel as there was not enough statistics to get a dependable 
ratio. 
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